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We compute go'D-n and gB'Sir using a framework in which all elements are constrained by Dyson- 
Schwinger equation studies of QCD, and therefore incorporates a consistent, direct and simultaneous 
description of light- and heavy-quarks and the states they may constitute. We link these couplings 
with the heavy-light-meson leptonic decay constants, and thereby obtain go'Dir ~ IS.Olj J and 
gB'Bn = 30.0l5;4. From the latter we infer gs ~ 0.371q;q2. A comparison between go'Dn and 
gs'Sw indicates that when the c-quark is a system's heaviest constituent, AqcD/jTic-corrections are 
not under good control. 
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Introduction. Non-perturbative QCD effects in 
heavy flavor physics arc a persistent challenge to an accu- 
rate determination of the Cabibbo-Kobayashi-Maskawa 
matrix elements, to virhich the BaBar, Belle, CDF, CLEO 
and FOCUS collaborations have dedicated significant ef- 
fort in the past decade. Besides improvements in the 
determination of these Standard Model parameters, the 
CLEO collaboration also performed a first measurement 
of the width of the D meson's nearest resonance. Their 
reported value for the charged vector meson, T{D*^) = 
96 ± 4 ± 22 keV [1], is of great interest because it opens 
a vifindow on nonperturbative strong physics involving 
heavy quarks. More specifically, it allows an extraction 
of the coupling go-'D-K, which some relate to a putative 
universal strong coupling, g, between heavy-light-vector 
and -pseudoscalar mesons and a low-momentum pion in a 
heavy-meson chiral Lagrangian [2]. The step from go'Drr 
to g is contentious, however, because the c-quark is not 
truly heavy and hence corrections to the heavy-quark 
limit are not necessarily under good control. 

In attempting to compute g, one may work with the 
matrix element 



{H{p2)iT{q)\H*{pi,X)) ^ gH'Hn ex-q 



(1) 



which defines the dimensionless coupling of a heavy-light 
vector meson, H*, characterized by a polarization state 
A, and a heavy-light pseudoscalar meson, H, to a soft 
pion with momentum q = pi — P2- This matrix element 
describes the physical processes D* — > Dtt, with both 
the final pseudoscalar mesons on-shell. It also serves 
in computing the unphysical soft-pion emission ampli- 
tude B* — > Bt: in the chiral limit, which defines gs'Sn- 
A theoretically consistent comparison between these two 
couplings can provide a quantitative indication of the de- 
gree to which notions of heavy-quark symmetry may be 
applied in the charm sector. 
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18 ±3 


0.61 ±0.10 


32 ±5 


0.40 ±0.06 



TABLE I. Calculated values of H* — >■ Hn couplings com- 
pared with experiment and other estimates. For the lattice- 
QCD results: the subscript indicates the number of dynam- 
ical light-fermions employed in the computation; the valence 
c-quark is treated directly but its dynamics is quenched in all 
simulations; and the B-meson simulations treat the heavy- 
quark as static. NB. Where useful, we have combined errors 
in quadrature in order to simplify presentation. Experimen- 
tally [3] (in GeV): m^o = 1.865, 771^.+ = 2.010, nigo = 5.280, 
ms- = 5.325, = 0.1396, f^+ = 0.1307. 

Selected results for gH'H-n and the associated value of 
gwH-K 
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(2) 



are listed in Table I. As we explain subsequently, no 
entirely ab initio approach to such hadronic decays is 
currently available. This might explain the convergence 
of modern theory for go* Dtt in the presence of an experi- 
mental result but the disagreement over gB'Bir, which is 
kincmatically forbidden. 
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Present-day simulations of lattice-regularized QCD 
treat the valence c-quark directly as a propagating mode 
but its dynamics is quenched, whereas the &-quark is con- 
sidered as static. Beyond statistical, there are errors 
associated with chiral extrapolation, discretization and 
perturbative renormalization of the currents involved. It 
is apparent from Table I that results obtained within this 
approach exhibit a significant difference between go and 
cjB- This is probably because c-quark physics is not well- 
approximated by the heavy-quark limit and hence related 
observables contain significant AqcD/mc corrections. 

QCD sum-rules, too, have been used to estimate 
gH'HTT- In advance of the CLEO measurement, there 
were numerous results, which differed amongst them- 
selves by as much as a factor of two; and all of which 
were more than 25% below the experimental value deter- 
mined subsequently [9]. The value obtained in this ap- 
proach depends heavily on the functional forms used for 
extrapolation to the mesons' on-shell momenta from the 
domains upon which the calculations are actually valid 
[10]. Based on Ref. [10], loop corrections in D-meson ef- 
fective field theory have been used to constrain extrapo- 
lation of the D*Dtt coupling to physical momenta [11]. 

The simultaneous computation of both couplings has 
been completed in a dispersive quark model [12]. How- 
ever, this approach makes no qualitative distinction be- 
tween light- and heavy-quarks. Hence it cannot ve- 
raciously describe dynamical chiral symmetry breaking 
(DCSB), which makes problematic its treatment of the 
structure and interactions of pions. 

Symmetry-preserving models built upon robust predic- 
tions of QCD's Dyson- Schwinger equations (DSEs) [IS- 
IS] provide a sound framework within which to examine 
heavy-meson observables [16-19]. Such studies describe 
quark propagation by fully dressed Schwinger functions. 
This has a material impact on light-quark characteristics 
[15]. Constrained by experimental and theoretical heavy- 
light-meson information available at the time, Ref. [18] 
obtained go'Dn = H, go = 0.37 and gB'Brc ~ 23, 
(jB = 0.29. These values compare well with a repre- 
sentative average of theoretical estimates then available; 
viz., go'DTT = 12 ± 4, gB'Bir = 25 ± 7 [20]. However, 
go'D-n = 11 is not consistent with the subsequent CLEO 
measurement [1] . Numerous advances have been made in 
experiment and theory in the intervening period. Herein, 
therefore, we reassess the DSE study [18], implement 
improvements detailed in Ref. [19] and add one more; 
namely, an expression of the probable difference in size 
between the vector and pseudoscalar heavy-light mesons. 

Framework and tools. At leading-order in a system- 
atic, nonperturbative, symmetry-preserving DSE trunca- 
tion scheme [21, 22], the H* — >■ Htt decay amplitude is 
given by 



gH'H,, £ ■ q = trcD 



(2^ 



• TH'{k;pi)SQ{kQ) 



X Tnik; -p2)Sfik'f)r^{k-, -g)5/(fc/) , (3) 
where: the trace is over color and spinor indices; kg — k-\~ 




k + Wipi - P2 



FIG. 1. (Color online) Pictorial representation of Eq. (3), our 
impulse approximation to the H* — >■ Hit decay: solid lines 
- dressed-quark propagators [Eqs. (9), (10), (12)]; and filled 
ellipses - meson Betlie-Salpeter amplitudes [Eqs. (13) - (15)]. 



wiPi, k'j: = k + wipi—p2, kf ~ k~W2Pi, and the relative- 
momentum partitioning parameters satisfy wi + W2 = 1; 

is the vector-meson polarization four-vector; and S 
and r, described below, are dressed-quark propagators 
and meson Bethe-Salpeter amplitudes, respectively. The 
approximation represented by the diagram in Fig. 1 has 
been widely used successfully; e.g., Refs. [13, 14, 23-25]. 
It is reasonable to expect that corrections associated with 
final-state interactions arc small owing to the large c- 
quark mass. 

Along with the radiative-pion couplings, we simulta- 
neously calculate the H*- and _ff-nicson Icptonic decay 
constants, which arc determined via [18, 26, 27]: 



/d*k 



7m X^.H' (fc; P) 



(4) 



(5) 



where x(fc; P) % {k + wiP)T{k; P)Sf,_ (k - W2P). The 
Bethe-Salpeter amplitudes are canonically normalized; 
i.e., 
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(7) 



with an analogous expression for the H* . 

For any quark flavor, the dressed-quark propagator has 
the general form 

S{p) = -ij-p<Jv{p^)+<Js{p^) = l/[ij-pA{p^) + B{p^)], 

(8) 

and can be obtained from QCD's gap equation [13]. In 
connection with light-quarks, it is a longstanding pre- 
diction of DSE studies that both the wave function 
renormalization, Z{p^) ~ 1/A{p'^), and dressed-quark 
mass, MijP') = B{p^) / AijP'), receive strong momentum- 
dependent corrections at infrared momenta: Z{p'^) is sup- 
pressed and Af(p^) enhanced. These features are an ex- 
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pression of DCSB and, plausibly, of confinement.^ The 
enhancement of M{p'^) is central to the appearance of a 
constituent-quark mass-scale and an existential prereq- 
uisite for Goldstone modes. These DSE predictions are 
confirmed in numerical simulations of lattice-QCD [28]. 

The impact of this infrared dressing on hadron phe- 
nomena has long been emphasized [23] and, while nu- 
merical solutions of the quark DSE are now readily ob- 
tained, the utility of an algebraic form for Siji) when 
calculations require the evaluation of multi-dimensional 
integrals is self-evident. An efficacious parametrization, 
which exhibits the features described above, has been 
used extensively; e.g., [18, 19, 29, 30]. It is expressed via 



+T{bix) T{b3x) [bo + 63 J'lea;)] 
1 



avix) 



[1- J'(2(a; + m2))] , 



(9) 
(10) 



with 



m/A, T{x) = [1 — exp(— x)]/x. 



tTs(x) = Xasip^) and av{x) = X^av{p^)- The param- 
eter values were fixed in Ref. [18] by requiring a least- 
squares fit to a wide range of light- and heavy-meson 
observables, and take the values:^ 



/ 



rrif 



d 0.00948 0.131 2.94 0.733 0.185 



(11) 



The mass-scale A = 0.566 GeV, with which value the 
current-quark mass is nid — 5.4 MeV. In addition one 
obtains the following Euclidean constituent-quark mass, 
defined as IVP = {s\s + M^{s) = 0}: Md = 0.36 GeV. 

As noted elsewhere [31, 32], studies which do not 
or cannot implement light-quark dressing in this QCD- 
consistent manner, invariably encounter problems be- 
cause of the need to employ rather large constituent- 
quark masses and the associated poles in the light-quark 
propagators. This translates into considerable model- 
sensitivity in the results for any heavy-light form factors, 
such as B- to light-meson transition form factors [33] . 

Whilst the impact of DCSB on light-quark propagators 
is significant, it is less so for heavier quarks. This is plain 
in Fig. 1 of Ref. [18]. It can also be seen by considering 
the renormalization-point-invariant ratio <;/ := af/Mj", 
where ct/ is a constituent-quark tr-term [14]. This ratio 
measures the effect of explicit chiral symmetry break- 
ing on the dressed- quark mass-function compared with 
the sum of the effects of explicit and dynamical chiral 
symmetry breaking. Calculation reveals [14] <;d = 0.02, 
<;s = 0.23, <jc = 0.65, c^b ~ 0.8, results which arc read- 
ily understood. Naturally, <;f vanishes in the chiral limit 



^ Equations (9), (10) represent S{p) as an entire function. This 
entails the absence of a Lohmann representation, which is a suf- 
ficient condition for confinement [15]. 

^ e = 10~* in Eq. (9) acts only to decouple the large- and 
intermediate-p'^ domains [23]. 



and must be small for light-quarks because the magni- 
tude of their constituent-mass owes primarily to DCSB. 
For heavy-quarks, <;f approaches unity because explicit 
chiral symmetry breaking becomes the dominant source 
of their mass. We therefore use a constituent-quark-like 
propagator for c- and fo-quarks; viz.. 



SQ{k) 



1 



ij-k + Mq 



Q 



(12) 



where [18]: A'h = 1.32 GeV, A/^ = 4.65 GeV. 

The meson Bcthe-Salpetcr amplitudes, which appear 
in Fig. 1 and are consistent with the generalized impulse 
approximation, are properly obtained from an improved- 
ladder Bethe-Salpeter equation [13]. The solution of this 
equation requires a simultaneous solution of the quark 
DSE. However, since we have already chosen to sim- 
plify the calculations by parametrizing S{jp)^ we follow 
Ref. [18] and also employ that expedient with F. The 
axial-vector Ward-Takahashi identity and DCSB have 
an enormous impact on the structure and properties 
of light pseudoscalar mesons. Indeed, the quark-level 
Goldberger-Treiman relations derived in Ref. [26] moti- 
vate and support the following efficacious parametriza- 
tion of the pion's Bethe-Salpeter amplitude: 



where := B„l.„_ 



i75^5x(fc2), 



(13) 



is obtained from Eqs. (8) - 
(10) through the replacement b"^ b^ = 0.204 [18]. 
Equation (13) expresses an intimate connection between 
the leading covariant in a pseudoscalar meson's Bethe- 
Salpeter amplitude and the scalar piece of the dressed- 
quark self-energy [26]. 

Whilst the renormalization-group-improved rainbow- 
ladder DSE kernel is appropriate for the study of mesons 
constituted from equal-mass constituents [13, 35], this is 
not so for heavy-light mesons. In such systems cancela- 
tions, which largely mask the effect of dressing the quark- 
gluon vertices, are blocked by the dressed-propagator 
asymmetry; e.g., a recent analysis [36] obtained a fair de- 
scription of D- and i?-meson masses but underestimated 
their leptonic decay constants by 30-50%. (The approach 
introduced in Ref. [37] might be the advance needed 
to make progress with heavy-light systems.) IMoreover, 
as we have already chosen to simplify the calculations 
by parametrizing Siji) and F^, it is rational to fol- 
low Rcfs. [18, 19] and employ that expedient, too, with 



exp(-fcVa;|f) 



Ml 



H 



^ _ exp (-fcVcj|f,) 



(14) 
(15) 



The normalization, A/'^(.), is fixed by Eqs. (6) and (7). 
Herein, however, we depart from Refs. [18, 19] and do 
not assume heavy-quark symmetry to be realized exactly; 
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namely, we eschew the spin-independent Ansatz and fit 
each of the parameters loh and ujh* to their respective 
leptonic decay constants, either known from experiment 
or estimated via lattice-simulations. [NB. Our results 
are not materially affected by the pointwise form of the 
functions in Eqs. (14), (15).] 

Poincare invariance is a feature of the direct applica- 
tion of DSEs to the calculation of hadron properties. 
However, that is compromised if one does not retain 
the complete structure of hadron bound-state amplitudes 
[27]. This applies herein because we use one-covariant 
models for the amplitudes. To proceed we must there- 
fore specify the relative momenta. 

As indicated in Fig. 1, in computing gH*H-n'- we al- 
locate a fraction wi of the heavy-light-meson's momen- 
tum to the heavy-quark; W2 to the light-quark; and 
momentum-conservation specifies the remaining momen- 
tum. An analogous procedure is followed for the lep- 
tonic decays and normalization, Eqs. (4) - (7). The 
relative-momentum partitioning parameter is defined via 
a center-of-mass prescription; viz.. 



Mh + Md 



(16) 



which yields = 0.79, — 0.93. In addition, we 
interpret the amplitudes T in Eqs. (13) - (15) as the 
zeroth Chebyshev moment of the leading-covariant in 
the relevant meson's Bethe-Salpeter amplitude. We em- 
phasize that any sensitivity to a definition of the rel- 
ative momenta is an artifact owing to our simplifica- 
tions. Every study that fails to retain the full struc- 
ture of the Bcthc-Salpctcr amplitude shares this compli- 
cation, which is never encountered in complete studies; 
e.g., Rcfs. [27, 35, 36]. 

Results. In Ref. [18] it was noted that it is a poor 
approximation to assume heavy-quark symmetry for c- 
quark systems because corrections can be as large as a 
factor of two in semileptonic c — d transitions. This is 
also seen in Ref. [19], wherein the assumption ujjj-^ = lud 
yields fo' - fo = 321 - 223 = 98MeV, whereas lattice- 
QCD combined with experiment suggests a smaller differ- 
ence: 39 ± 29 MeV. In relation to 5-quark mesons, whilst 
lob* = is more likely to be a good approximation, 
it cannot be exact and hence it is worth exploring the 
impact of relaxing this constraint. 

Wc therefore have four parameters: lod, i-od*, and 
a;^., which we determine via a least-squares fit to the 
leptonic decay constants [in MeV]: /n = 205.8 ±8.9 [38]; 
fo' = 245 ±20 [39]; /b = 229^^^ [40]; Jb- = 1961^ 
[39], using Eqs. (4) and (5). A perfect fit is possible, and 
is obtained with [in GcV]: 



ujb 



1.26 ±0.10. 
2.26 



+0.76 
-0.68 ■ 



UJD' 
UJB- 



1.21 ±0.14, 
1.58 



hO.52 
-0.27 ■ 



(17) 



It is notable that, in contrast to Ref. [19], the 
width parameters are all consistent with intuition: 
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FIG. 2. (Color online) Dimensionless coupling go'Div'- com- 
parison between experiment and recent estimates, with their 
associated errors added in quadrature. (See Legend and Ta- 
ble I for details.) 



:= 1/cjB = 0.09fm < = 0.12 fm < = 

0.157 fm < £u]*^ = 0.163 fm; namely, by this rudimen- 
tary gauge, pseudoscalar mesons are smaller than vec- 
tor mesons and systems containing a single 6-quark are 
smaller than those containing a c-quark. 

Using the width parameters in Eq. (17), we computed 
the strong i7*+i/°7r+ couplings 

go'D^ = fim go'Diriq^) , gB'Bi, ^ fim gB'Buiq^) , 

(18) 

where the latter serves as a definition since the process 
B* — > B-K is kinematically forbidden. We list our results 
in Table I and, for go'Din present a pictorial comparison 
with experiment and other studies in Fig. 2. Importantly, 
in our approach one can directly calculate the amplitude 
at the fully on-shell point and with the physical light- 
quark current-mass: no extrapolations are necessary. 

There are two obvious sources of uncertainty in our re- 
sults for gH'HTT- The errors in the values of the leptonic 
decay constants translate into the uncertainties quoted 
on the width parameters in Eq. (17) and produce a range 
of values for gH-H-K- In addition, a change of ±20% in 
gives variations in go'D-K and gB*B-K of ~ 10% and 



c,b 



~ 5%, respectively. We treat these variations as inde- 
pendent uncertainties on each coupling and add them in 
quadrature to produce the errors quoted in Table I. 

Our results are significantly different to those re- 
ported in Ref. [18] because, as noted above, there are two 
material differences between this calculation and that. 
Namely, we allow: lod* ^ lod and lob' ^b] and the 
light-quark to carry a fraction of the i/*-meson's mo- 
mentum. In connection with gB'B-m since the 6-quark is 
genuinely heavy and should therefore carry most of the 
S*-meson's momentum, the improvement arises primar- 
ily because ujb' ^b- On the other hand, the c-quark 
is neither light nor truly heavy and thus go-DTi is quite 
sensitive to the amount of the D*-meson's momentum 
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carried by the light-quark. If we require W2 = 0, as in 
Ref. [18], then one finds gD*D-n < 13; even allowing for 
ijJy)- ^ lotj. a less important factor is our use of updated 
D, £>*, B and B* masses. These values influence the nor- 
malization constants, Eq. (7), but not to an extent which 
requires the effect of experimental mass uncertainties to 
be included in our error estimate for gn-'HTr- 

Summary and Conclusions. We presented a cal- 
culation of qd'Ott and gs-'BTr based upon QCD's Dyson- 
Schwinger equations (DSEs). By implementing a more 
realistic representation of heavy-light mesons; e.g., al- 
lowing the light-quark in the heavy-light-meson to carry 
some of the meson's momentum and for a difference be- 
tween the sizes of pseudoscalar and vector mesons, our 
analysis improves significantly upon earlier DSE-based 
studies. Furthermore, we step beyond other models be- 
cause, by expressing confinement and dynamical chiral 
symmetry breaking (DCSB) in a manner compatible with 
the predictions of QCD's DSEs, our approach incorpo- 
rates a consistent, direct and simultaneous description 
of light- and heavy-quarks and the states they may con- 
stitute. Finally, even with respect to modern numerical 
simulations of lattice-QCD, our approach has merits; e.g.. 
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